Importance of subleading corrections for the Mott critical point 
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The interaction-induced metal-insulator transition should be in the Ising universality class. Ex- 
periments on layered organic superconductors suggest that the observed critical endpoint of the 
first-order Mott transition belongs instead to a different universality class. To address this ques- 
tion, we use dynamical mean-field theory and a cluster generalization that is necessary to account 
for short-range spatial correlations in two dimensions. Such calculations can give information on 
crossover effects, in particular quantum ones, that are not included in the simplest mean-field. In the 
cluster calculation, a canonical transformation that minimizes the sign problem in continuous-time 
quantum Monte Carlo calculations allows us to obtain very accurate results for double occupancy. 
These results show that there are important subleading corrections that can lead to apparent expo- 
nents that are different from mean-field. Experiments on optical lattices could verify our predictions. 



Half-filled band materials should be metallic, but they 
are sometimes insulators Ij. This paradox was discussed 
by Boer and Verwey and by Peierls as early as 1937, but 
the first theoretical advance came from Mott in 1949. 
He found that as a function of some external parameter, 
it is possible to control the ratio of interaction energy 
to kinetic energy and drive the system through a metal- 
insulator transition. This Mott transition has by now 
been clearly identified in a few materials [1 and in optical 
lattices of cold atoms [31 The order parameter for 
the interaction-induced transition should be in the Ising 
universality class 4-6 , with no breaking of translational 
or rotational invariance. This has been verified explicitly 
in the three-dimensional compound V2O3 [7]. 

It thus came as a surprise when it was discovered that 
in two-dimensional layered k-ET organic superconduc- 
tors [8], critical exponents for the Mott critical point, 
measured in both charge (conductivity) [S] and spin 
(NMR) channels [TU], appear to belong to a different 
universality class. Several proposals have appeared to 
explain this result. Imada et al. [TTl [12] suggested that 
while the high-temperature regime is described by classi- 
cal Ising exponents, there is also a continuous transition 
at T = and, in between, a marginal quantum critical 
point that controls the observed behavior. Papanikolaou 
et al. [13] instead started from the 2d-Ising universal- 
ity class and argued that, away from criticality, the sub- 
leading energy exponent dominates for the conductivity 
over the leading order parameter exponent. The latter 
becomes relevant only very close to Tc. A recent experi- 
ment on thermal expansion coefficient finally, argues that 
the Ising universality class is the correct one [14] . That 
finding disagrees with the latest theoretical calculation 
[T5] performed with Cluster Dynamical Mean-Field the- 
ory (CDMFT) [m [TT] that measured an exponent 5 = 2 
in agreement with the above-mentionned conductivity [5] 
and NMR experiments [lO] . 

Here we revisit the critical behavior at the Mott critical 
endpoint by studying the one-band Hubbard model, the 



simplest model of interacting electrons that contains the 
physics of the Mott transition. We use single-site Dynam- 
ical Mean-Field Theory (DMFT) and CDMFT. 
Single-site DMFT is exact in infinite dimension and can 
be applied to lower-dimensional lattices '51] [52] . While it 
can be proven analytically j5j that for DMFT the behav- 
ior is mean-field like, the possibility of quantum critical 
transients, the size of the critical region and the precise 
value of the exponents have not been verified numeri- 
cally. This is important since there could be an unstable 
quantum critical point controlling the behavior in the ex- 
perimentally accessible regime, as has already been ob- 
served for the conductivity |23| . The same question arises 
with CDMFT that takes into account the momentum de- 
pendence of the self-energy, a physical ingredient that is 
known to be important in two dimensions |24H35j . Since 
the Mott transition occurs far above long-ranged ordered 
ground states, CDMFT should be an accurate description 
of the organics, except asymptotically close to the tran- 
sition where the usual critical fluctuations should take 
over mean-field behavior. 

Although this problem has already been studied with 
these methods, improvements in computer performance 
and in algorithms allow us to obtain much more accu- 
rate data. In the case of CDMFT, for the frustrated 
lattice considered here, the sign problem in the Contin- 
uous Time Quantum-Monte Carlo solution of the Hy- 
bridization expansion (CT-HYB) |36H39j is minimized by 
a canonical transformation. This allows us to approach 
the critical point ten times closer in reduced pressure 
than previously possible. 

Method: The simplest model that contains both the 
strong on-site Coulomb repulsion and the kinetic energy 
of the frustrated k-ET's lattice, is the half-filled Hubbard 
model on a 2D anisotropic triangular lattice 
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where creates a spin a electron at site z, n^o- = c'ia^ia- 
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FIG. 1: Periodic partitioning of the anisotropic triangular 
lattice into 2x2 plaquettes for CDMFT. 

is the spin a density at site i, tij = t*^ are the hopping 
amplitudes as shown in Fig. ([I]) while /i and U are, respec- 
tively, the chemical potential and the screened Coulomb 
repulsion. 

We use single-site DMFT j21j and its cluster extension 
CDMFT [iniin] to solve the Hamiltonian Eq.Q. These 
methods start with a periodic partitioning of the infinite 
lattice model into independent sites (DMFT) or clusters 
(CDMFT). The missing environment of the cluster is re- 
placed by a bath of non-interacting electrons. The action 
of the cluster in a bath model may be written as 

S' = 5ci(c1',c)+ / dTdT'c^T')A{T' -t)c{t), (2) 

where Sd is the cluster-action as obtained by the parti- 
tioning, c the column vector of the corresponding Cio-'s 
and the bath has been integrated out in favor of a hy- 
bridization function A — (Aia-j^i ) . This defines an effec- 
tive impurity model. Approximating the unknown lattice 
self-energy locally by the impurity self-energy S' and us- 
ing Dyson's equation, the lattice Green's function reads 

with Gojat the non-interacting lattice Green's func- 
tion. The hybridization function A is determined self- 
consistently from the requirement that the impurity 
Green's function computed from the action Eq. ^ co- 
incides with the projection of the approximate lattice 
Green's function G[^^^ on the impurity cluster. This self- 
consistency condition may be derived from a variational 
principle gD]. For CDMFT we take the 2 x 2 plaquette 
illustrated in Fig. ([!]). This accounts at least locally for 
the geometrical frustration in the k-ET. 

To obtain the impurity Green's function (and other 
observables), we use a continuous time quantum Monte- 
Carlo (CTQMC) solver [SnHSS]- This method starts with 
a diagrammatic expansion of the impurity partition func- 
tion in powers of the hybridization function. The key idea 
is then to recombine, for a given 'cluster-vertex', all pos- 
sible contractions over the bath (i.e. the hybridization 
function) by a determinant before summing up the series 
by Monte-Carlo sampling. Without this recombination 



of diagrams the fermionic sign problem would be fatal 

m- 

In the case of CDMFT, symmetries of the problem 
can be used to speed up the simulation by choosing a 
single particle basis in Eq. ^ that transforms according 
to the irreducible representations [39 . In our case, sepa- 
rate charge conservation of a =tj i particles and the C2v 
point group symmetry of the anisotropic plaquette lead 
to the single particle basis (see Fig. [l] for indices) 

CBia = ~ '^3<j) (4) 

with Ai,Bi and B2 irreducible representations of C2V 
{A2 is empty). Due to the degeneracy in the Ai subspace, 
there is a degree of freedom in the choice of basis which 
may be parameterized by an angle 9 as follows, 

cos9c'j^^„ - s'mOcAicr, s'mOc'^^^ + cos Oca,^^- (5) 

In this basis the hybridization function A takes a block- 
diagonal form with one 2x2 block (Ai) and two 1x1 
blocks {Bi and B2) for each spin (in the normal phase). 
The sign problem in the Monte Carlo simulation shows a 
strong dependence in as shown in Fig. [2] for t/t' = 0.8, 
/3 = 20 and different values of U. One can check that the 
maximum of the average sign is related to the angle 9 that 
minimizes the off-diagonal elements of the hybridization 
function (^1 block) with respect to some norm. The dots 
in the inset of Fig. [2] indicate the maximum with respect 
to Li and L2 on [0, /?]. The usual basis, 9 — 0, has a bad 
sign problem. 

To analyze our results, we derive subleading correc- 
tions to mean-field theory. The singular part of the 
mean-field equation for the order parameter rj takes the 
form [Smi] 

prj + af — h (6) 

with c a constant, while p and h are defined by 
p = Pi{U- Uc) +P2{T- Tc) and h = hi{U- U^) + 
h2 (T — Tc) ■ Like in the liquid-gas transition, interaction 
strength and temperature are not in general eigendirec- 
tions, which explains the way they appear in p and h. 
When p = 0, the solution is 77 = {h/cY^^ , which defines 
6 — 3. Approaching the critical line along 6U = [U — Uc) 
for example, the mean-field Eq.([6]) takes the form 

Pi6Ur] + CT]^ = hiSU. (7) 

One can show that the general solution of that equation 
is of the form 

00 

r, = Y,SU'^\. (8) 

i=l 
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FIG. 2: Average sign in CTQMC simulations of the 
anisotropic plaquette impurity problem at CDMFT self- 
consistency with t/t' — 0.8 {t = 1) and /3 = 20 as a function 
of the angle 9 in Eq. [5] for different values of U. The inset 
zooms on the region where the sign takes its maximum, as 
indicated. The dots associated with each curve indicate the 
angle where the off-diagonal elements of the corresponding 
hybridization functions are minimal with respect to the Li 
norm (filled) and the L2 norm (empty). 

with expansion coefRcients rji. The first term, SU^^^, and 
the subleading correction, SU^^^, are the only terms that 
lead to an infinite first derivative at the critical point. 
In the case of DMFT, 77 is the singular part of the hy- 
bridization function. 

In the following we compute double occupancy D = 
(n^ni) . Double occupancy in general should be a smooth 
function of rj that can be expanded as a power series, a 
result that can be proven in DMFT [5\ . Hence, even when 
77 is dominated by the leading term SU^^^ , the rj"^ term of 
the power series leads to subleading 5U^^^ corrections. 

For the single-band Hubbard model, singular be- 
havior of D implies singular behavior in both spin 
and charge channels [5], as follows from the follow- 
ing two sum rules on spin, Xsp, and charge, Xchi sus- 
ceptibilities, TJ2nI ■i§^Xsp{<l,^n) = n - 2D and 

/ -(^Xch (qj'^n) = n + 2D-n'^ where w„ are Mat- 
subara frequencies and q wave vectors in the Brillouin 
zone. 

Below the critical temperature, there is a first-order 
transition with a jump in double occupancy that scales 
like with /3 — 1/2. It is very difficult to obtain this 
exponent numerically because of hysteresis. Similarly, 
the exponent for the susceptibility {dij/dh)^ ^ p'"' with 
7 = 1 requires numerical differentiation and cannot be 
obtained accurately. 

Results: Figure [3] displays double occupancy as a func- 
tion of interaction strength calculated for both single-site 
DMFT (blue squares) and CDMFT (red circles) at our 
best estimate of the corresponding critical temperatures. 
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FIG. 3: Double occupancy as a function of U near the Mott 
critical point for the Hubbard model on an anisotropic tri- 
angular lattice with t' /t = 0.8 {t = 1) at half-filling and 
fixed critical inverse temperature j3 = 11.15 (squares) for 
DMFT and /3 = 9.9 (circles) for CDMFT on a 2 x 2 plaque- 
tte. The sohd fines show a fit to f{U) = cisgn{SU)\5U\^^^ + 
C2\5U\'^/'^ + C3SU + Dc {5U = U - Uc) with the same pa- 
rameters ci, C2, C3, Dc,Uc and 5 for the metallic (filled sym- 
bols) and the insulating region (open symbols). The best 
fitting values {Uc,Dc,d) are (10.445,0.0325,2.93) for DMFT 
and (7.9316, 0.0679, 3.04) for CDMFT. (a) Linear plot cen- 
tered at {Uc,Dc). The insets zooms on the regions close to 
the critical point, (b) Logarithmic plot in reduced units rela- 
tive to the critical point with CDMFT data shifted by a factor 
of 10^'^ along the y-axis. The dashed lines show the function 
(X \U — Uc\^''^ with 5 as indicated. In the critical regime, up 
to 500 iterations are necessary for convergence in the itera- 
tive solution of the (C)DMFT equation. Once convergence 
is reached, we take the average over hundreds of iterations. 
Monte-Carlo sweeps per iteration: 6 • 10^ for DMFT and 10^ 
for CDMFT. 

Both the metallic (filled symbols) and insulating (open 
symbols) sides are shown. The critical temperature is 
found as follows. Below the critical temperature, there 
is hysteresis and a jump in double occupancy. Above 
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the critical temperature, double occupancy is continuous. 
First we searched for the highest (lowest) temperature 
where hysteresis (continuity) can be checked in a reason- 
able time. The mean of these two temperatures is then 
taken as an approximation for the critical temperature. 
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0.32 


3.05 


0.33 


ImAioc(t^n = vr//?) 


3.02 
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TABLE I: Estimates of the exponent S from a fit of Eq. [9] 
to the critical behavior of the double occupancy D, the local 
Greens function Gioc at t = 13/2 and the real and imaginary 
parts of the local hybridization Aioc function at the lowest 
Matsubara frequency, as obtained by DMFT and CDMFT 
for the same model and parameters as in Fig. [3] The ratio 
|c2/ci| indicates the weight of the subleading correction, as 
seen from Eq [9] 

In Fig. the scale is linear whereas it is logarithmic 
in Fig. ^p. The solid lines are fits to the functional form 
suggested by Eq. ^ and by the smoothness hypothesis 
for D, 

D-Da = ciSgn{5U)\SU\^/^ + C2|(5C/|^/* + C3SU (9) 

where S and the coefficients are adjustable parameters. 
The fits include both the metallic and the insulating 
sides. We find S = 2.93 ± 0.15 for DMFT, where we 
know that the analytical result [5] is (5 = 3. For CDMFT 
we find S = 3.04±0.25. The errors are estimated from the 
values of d at the two temperatures just below and above 
the critical one. The log-log plot in Fig. [3)3 shows that 
the data does not lie on a perfect straight line over the 
wide range of reduced units considered here. The straight 
dashed lines are guides to the eye that show that the ex- 
ponent that we would obtain by fitting over a limited 
range of 6U would decrease from 6 = 3 towards 6 = 3/2 
as we move away from the critical point. On the metal- 
lic side, the crossover extends over a rather wide region 
where the the exponent is close to 6 = 2. 

As shown in Table |T] different critical quantities lead 
to coherent estimates of i5, whereas the importance of 
the subleading corrections varies strongly from case to 
case. Note that an asymmetry of the critical quantity 
between the metallic and the insulating sides, as visible 
in the double occupancy, indicates important subleading 
corrections. 

Discussion: We have shown that the results of nu- 
merical calculations with both DMFT and CDMFT are 
consistent with 6 = 3 (as predicted analytically [S] for 
DMFT), if we include the subleading corrections in the 
analysis. These are particularly important if the critical 



data is asymmetric in the accessible region. Otherwise, 
fits with a single exponent may lead to 5 closer to i5 = 2 
in the metallic phase, as observed in the organics [H [10] , 
and in previous CDMFT [TS] calculations, where ana- 
lytical results are not available. Note that in the latter 
calculations, the value of t' /t and the direction of ap- 
proach to the critical point differ from ours. That may 
modify the size of the crossover region. 

Extracting the pressure dependence of model param- 
eters from band structure calculations [H], we estimate 
that our numerical results are as close to the critical point 
in reduced units as are the experiments. The value 7 1 
in these experiments is the same as the mean-field one, 
while (3 = 1 would imply that a non-singular term domi- 
nates the physics in the accessible range. Since coupling 
to the lattice would favor mean-field exponents [33], it 
would be interesting to reanalyze the experimental re- 
sults by including the subleading correction to the mean- 
field behavior. If the coupling to the lattice is irrelevant, 
the exponents should eventually differ from mean-field 
behavior sufficiently close to the transition. 

To definitely settle this issue experimentally, it would 
be interesting to study the two-dimensional Mott tran- 
sition in frustrated optical lattices, where double occu- 
pancy is directly accessible [5]. Our derivation of the 
subleading correction also shows that it can possibly be 
relevant in more mundane cases when mean-field expo- 
nents are appropriate [44], away from the critical regime. 
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